Electron transport in interacting hybrid mesoscopic systems 
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A unified theory for the current through a mesoscopic region of interacting electrons connected to 
two leads which can be either ferromagnet or superconductor is presented, yielding Meir-Wingreen- 
type formulas when applied to specific circumstances. In such a formulation, the requirement of 
gauge invariance is satisfied automatically. Moreover, one can judge unambiguously what quantities 
can be measured in the transport experiment. 
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Mesoscopic electron transport has received an increas- 
ing attention both theoretically and experimentally in 
last decade.!]. In mesoscopic or nanoscale systems the 
wave nature of electrons becomes apparent and the 
transport process is coherent. The Landauer-Buttiker 
formula 2], which encodes the current in the local prop- 
erties of the interacting mesoscopic region and the equi- 
librium distribution functions of the noninteracting elec- 
tron reservoirs, enhances our understanding of meso- 
scopic electron transport and has been applied successful 
in many fieldst^]. In 1992 Meir and Wingreen 0] pre- 
sented a formulation for electron transporting through 
a small confined region (quantum dot, QD) where the 
electron-electron interaction is important, and recovered 
the form of Landauer-Biittiker formula in the noninter- 
acting case. 

Recent advances in nanofabrication and material 
growth technologies make it possible to realize various 
kinds of hybrid mesoscopic structures 0, H, S IS ^i 
nil of which the building blocks are normal met- 
als (N), ferromagnets (F) and superconductors (S). It 
is known that transport in the presence of a ferromag- 
net and a superconductor will be strongly related to the 
spin polarization of the ferromagnet and the Andreev re- 
flection at the boundary of the superconductor (isl ll^ . 
The co-existence of two ferromagnets or two supercon- 
ductors is revealed to display spin-valve effect 0| or 
Josephson effect H^. When an interacting normal metal 
is connected with bulk ferromagnet (s) and/or supercon- 
ductor (s), it is expected that the interplay among the 
electron-electron interaction, spin imbalance, Andreev 
reflection would induce more interesting or even more 
surprising features in mesoscopic electronic transport. 
Previous theoretical investigations on the transport prop- 
erties of specific mesoscopic hybrid structures, such as 
N-QD-S[ia, S-QD-SHi, F-QD-F[li, F-QD-S ^ etc., 
make some intuitive presumptions and hence lack math- 
ematical rigidity, which will be discussed below in detail 

In this paper, we provide a scheme to treat the trans- 
port problem in an interacting hybrid mesoscopic struc- 



ture in a unified way by using the Keldysh formalism |j, 
EH - It is shown that gauge invariance can be satisfied 
automatically. It is also shown that what physical quan- 
tities can be measured in experiment. 
We start with the Hamiltonian 



H = He + H-ji -I- He + Hr, 



(1) 



}^ is the 



Hamiltonian for the 



where Tic 

central interaction region. He (Hy ) fo r the left (right) 
lead can be either the Stoner model [l5j characterized by 
an exchange magnetization h with polar angle 9 or the 
5(75" Hamiltonian |_16J with order parameter A = jAje*"^: 
^ J^kA^fka + sign{a)h^ COS 9^ f - ^i',]fl,,„f^ka + 
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tunneling Hamiltonian 
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can be written 
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H.c. 



as 
or 



^kn\c 



jka'^dna- 

fJ-f/c is the chemical potential of the corresponding part, 
and hereafter the notations ct =ti i and cr = ± are used 
interchangeably. 

To see the tunneling processes more clearly, and even 
more importantly, to facilitate the analysis of the gauge 
invariance, and the simplification of the general current 
formula Q to the form of specific systems, we first trans- 
form the Stoner and BCS Hamiltonian and the tunnel- 
ing Hamiltonian by the following Bogoliubov transforma- 
tions 



H.c. 



Here -f ^ C,TZ, 



fyka = cos{9yf/2)ipyfka- - sgn{a) s\tl(9^ f /2)il) 



ifkcr 



sk'i^fsk'a + sgn{(j) sin 9^s'Pi''^ 



where 6. 



■ysk' 



arctan[(eT,fc„ + Je^^.^ + A'^)/{e^k 



A2)]i/2 and V (V^) is the pair destruction (cre- 
ation) operator guaranteeing the particle conservation. 
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transforming a state of a given N particles into that with 
(7V+2)/(7V-2) particles, i.e., VyV\N) = |7V+2)/|iV-2) 
, thus make sense the abnormal off-diagonal Green's func- 
tions consists of two creation or destruction particle op- 
erators. One finds that the lead Hamiltonian is diagonal- 
ized after the above Bogoliubov transformations. 

In order to treat the ferromagnet and su- 
perconductor on the same footing, we intro- 
duce a 4-dimensional Nambu-spinor space, de- 
noted by = {il^l^ Vai V'lx V'at)^ and 
the Green's function in the Keldysh formal- 
ism H Go^.pihM) - i(Tc(*a(tl) ® *lj(i2))), 
where Tc is the contour-order operator, including 

gS(*i.^2) = T*t?(±ii Ti2)({*a(ii),*^(i2)}+), and 

^I'jitlM - ±*(*;j(t2)/*a(il) ® *a(tl)/*^(t2)). 

The tunneling Hamiltonian in such a representation 
takes the form 



7(F) 
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T.i^'^\fkmmcln+H.c), (2) 

E(*U'V^:(t)*,„+i/.c.), (3) 



where WUlt) = r/( Y)V^,{P(a*,cO, ^Zit) = 
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and the unitary rotation and phase operators(matrices) 
are given by 



cosx 





sin a; 


N 





cosx 





— sinx 


— sin a; 





cosx 








sinx 





cosx / 


cos X 


— "Psinx 





V* sin a; 


cos X 













cosx 










~V* sinx 












e 


— ix/h 






























P(x) 



The chemical potential /i^/c is incorporated into the 
phase matrix P and the ferromagnetism and supercon- 
ductivity are reflected in the corresponding rotation ma- 
trices R-'' and Rf . The tunneling Hamiltonian now rep- 
resents the explicit physical processes in the semiconduc- 
tor model 16]: an electron of spin a in the central regime 

can tunnel into either the spin a band or a band of the 
ferromagnetic lead, or tunnel into a spin a state or con- 
densate into an electron pair with a hole state of opposite 



spin being created; and vice versa. The total probabil- 
ity of the two tunneling processes into the same lead is 
cos^ X -1- sin^ X = 1 . As shown below, these rotation and 
phase matrices are very useful in our analysis of gauge in- 
variance, and even more importantly, the simplification 
of the formulas. 

The current from the left lead into the interacting re- 
gion is'S] 



lc{t) = -e{Nc) 



2e. 



nk 



+ ^ts,s,nm{tM)G%^Mn{il,t)) , (4) 



where the self-energy matrices after converting the sum 
iiito an integral / (iefcp^y^''(efc) (where pj/^^^ is the 
spin-dependent /normal density of states of the feromag- 
net /superconductor) are 
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and [x = g, g* , p) 

/ x'^''{e-c) 
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x'^'^ (e — c) 
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Here /(a;) = 1/(1 + /ksT) is the Fermi distribution 
function. The quasi-particle density of states of the BCS 
superconductor is p^^{£) — Re{g^^{e)}. 

The current flowing from the right lead 2k can be ob- 
tained in a similar way. In the steady transport problem, 
the current is uniform if no charge piles up in the central 
region, that is, Ic =^ ^^n- Symmetrizing Eq. |0J one 
finds 



i=l,3 



— OO 



(9) 



where the trace is over the level indices in the central 
region. Eq. lO with the self-energy matrices H5lt)l7l8|) 
is the central result of this work, it expresses the cur- 
rent through an interacting region in terms of the local 
properties of such region (G^d^^) ^'^'^ ^^"^ equilibrium 
distribution functions of the attached leads, as 

in the work of Meir and Wingreen Q • It can be applied to 
many types of hybrid mesoscopic structures even in the 
non-equilibrium situation, allowing various kinds of in- 
teractions in the central region. Notice that, the current 
is generally time-independent, except the case of two su- 
perconducting leads with nonzero bias, and that the full 
Green's functions G^^^ should be evaluated with the 
consideration of the tunnelling between the interacting 
region and the leads. The retarded/advanced Green's 
functions G^^ can be calculated in several approaches, 
such as the equation of motion formalisr nl22l , interpola- 
tive method[23j and the NCA technique [2J|. While the 
lesser one G^^ can be obtained from the Keldysh's equa- 
tion based on Ng's Ansatz|2^. With the help of the uni- 
tary property of the phase operator P, it is not difficult 
to find that Eq. ^ is gauge invariant, namely, X{t) re- 
mains unchanged under a global energy shift. When the 



bias V = {^c ~ becomes zero, one can readily find 

that I{t) = except the case of two superconducting 
leads with different superconducting phases, due to the 
coherent transport of quasi-particle pairs |2^. 

In the following, we apply this generalized result Eq. 
© to the specific structures we are interested in. When 
the two leads are both ferromagnetic, after a phase and 
a rotation transformation, Eq. jnj reduces to 



de , 



rr|([f^/(eTey)-r^^(e) 



[^^^(eTey)f£(£Tel/)- 
^^/(e)fK(e)] {G^Ae) - G2,rf(£)]) ..},(10) 



where t^f = R/t(^)r^/R/(^), 



9£/ - 6*7?,/, and 
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J d(t-t')e-(*-*')/''P(M7^cOR^(^) 



One sees that Eq. (|1(J|I is formally the same as the Meir- 
Wingree formula |j| in the normal lead case. The cur- 
rent is time-independent as one might expect, and just 
depends on the relative angle 9 / between two magnetiza- 
tion orientations and the bias V, the difference between 
the chemical potential of the two leads. 

The current through an interaction region with a fer- 
romagnetic and a superconducting lead can be derived 
similarly 



. i=l,3 

< 



grr{([r^/(eTey)-rf (e) 

[T''f'ieTeV){c{eTeV)- 



, {e){nie)\[G,Ae)-Gl,ie)]) },(11) 



in which the full Green's functions are 

r,a/< 



G 



d.d 
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d{t-t')^-^'-''y^v{^inct + ^) 



RH'-f)G:i<MRf\'-f)P\,nct' + 

The current through a generic hybrid mesoscopic struc- 
ture also depends on just the bias V. The depen- 
dence on the magnetization orientation of the ferromag- 
netic lead and the phase of the order parameter of the 
superconductor lead is absent after the rotation and 
phase operations, which are just presumptions in pre- 
vious investigations [itIIiI Im l20| . 

When the two leads are both superconductors, the sit- 
uation becomes complicated. However, we still obtain an 
elegant formula in this case following the same procedure 
as in the above derivation 



j=l,3 

- F 
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fc{eTeV)-Tf%e)in{e)]Gl,{e;t)) }, (12) 



where 

'■^Cs 



pt(eVt + ^)T%{sTeV) 



with ifs = ^pc — ^v,- Here we have added the time vari- 
able t into the full Green's function G2^^^(£;t), other 

than G2^^^(e). The reason is that the full Green's func- 
tions should be calculated in the presence of tunneling as 
well as interactions in the central region. In the present 
case, the ^-dependence can not be avoided in the self- 
energy matrices (I6I8II . while it can be removed by a uni- 
tary phase operation in the case of only one superconduc- 
tor lead. The current through an interacting mesoscopic 
region with two superconductor leads, is generally time 
dependent, as in the case of weak Josephson links [lit. 
However, in the limiting case of zero bias, the current 
is a time-independent nonzero quantity, as can be found 
fromEq. ifT^ . 

The time-dependence of G^^^^ arises from the cou- 
plings to the voltage biased two superconducting leads. 
One can show that G^'^^*^ depends only on the single 
time variable t at least within perturbative analysis. In 



fact the Green's function G 



r,a/< 



can be expanded in 



powers of the fundamental frequency lvq = 2eV/h, i.e., 
- E™G2;/<(e,e + mu;o/2)e*™"«*/2, which 

with the expression for the Green's function G^^^ be- 
low Ea. ((T^ is exactly the form of the double-energy 
transformation[2^ . And the current can be generally 
expressed as J{t) = Yl,n -^ne*""°*. 

The ferromagnetic or superconductor lead will be in 
a normal state, when the magnetization h or the or- 
der parameter A becomes zero. Such an observation 
allows us to study a more broad category of the hy- 
brid structures using the above formalism. By defining 
the unitary operators, we prove rigorously that the rel- 
ative value between some quantities, such as chemical 
potential, magnetization orientation, and order parame- 
ter phase etc. can be measured. All these quantities can 
he expressed as the energy-independent arguments of the 
exponential Junctions or the triangle functions in a uni- 
tary matrix^ reflecting some kind of requirement of the 
symmetrical invariance. As a further example, we re- 
cover all the known formulas obtained in the absence of 
interactions^ El El 113, El . This proves the validity 
and generality of our formalism in other respects. Taking 
into consideration the electron-electron interactions, one 



expect rich physics to show up, especially in the Kondo 
regime. 

In summary, we have given a unified formula for the 
current through an interacting region with either super- 
conducting or ferromagnetic leads. The current formula 
derived satisfies gauge invariance automatically. Such a 
current formula can be applied to an appreciable class 
of hybrid mesoscopic systems, allowing arbitrary interac- 
tions within the central nanoscale region. 
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